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Using optical-pump terahertz-probe spectroscopy, we study the relaxation dynamics of photoex- 
cited carriers in graphene at different temperatures. We find that at lower temperatures the tail 
of the relaxation transients as measured by the differential probe transmission becomes slower, ex- 
tending beyond several hundred picoseconds at temperatures below 50K. We interpret the observed 
relaxation transients as resulting from the cooling of the photoexcited carriers via phonon emission. 
The slow cooling of the photoexcited carriers at low temperatures is attributed to the bulk of the 
electron and hole energy distributions moving close enough to the Dirac point that both intraband 
and interband scattering of carriers via optical phonon emission becomes inefficient for removing 
heat from the carriers. Our model, which includes intraband carrier scattering and interband carrier 
recombination and generation, agrees very well with the experimental observations. 



In recent years, graphene has gathered much interest 
for electrical and optical applications due to its unusual 
band structure and properties^. Graphene has been in- 
volved in an array of applications that highlight its ver- 
satility and novelty as a platform for electronic, plas- 
monic, optical/IR, and terahertz devices^^. Realization 
of many of the graphene based devices relies on a good 
understanding of the nonequilibrium carrier and phonon 
processes and their associated time scales. In particular, 
the dynamics associated with the cooling and recombi- 
nation of photoexcited carriers are of interest in demon- 
strated and proposed graphene optoelectronic and tera- 
hertz devices^H^. 

Relaxation dynamics of photoexcited carriers in 
graphene have recently been studied using ultrafast opti- 
cal/IR pump-probe spec troscopy by several groups in- 
cluding the author a 9 " 16 !. These measurements, which 
typically measure the relaxation of the high energy tail 
of the carrier distribution, have shown that the photoex- 
cited carriers thermalize within few tens of femtoseconds 
to generate a hot carrier distribution. This hot distribu- 
tion then cools rapidly via optical phonon emission on 
a time scale of hundreds of femtoseconds. Within one 
picosecond, the carrier and the optical phonon tempera- 
tures equilibrate, and carrier cooling slows. At this point, 
cooling is limited by the exchange of energy between the 
carriers and the optical phonons and the subsequent an- 
harmonic decay of optical phonons into acoustic phonons. 
The questions that still remain unanswered pertain to the 
nature of the relaxation dynamics over much longer time 
scales and to the role played by carrier generation and 
recombination processes in the observed relaxation dy- 
namics. The answer to these questions is interesting both 
from the perspective of practical devices and also from a 
theoretical point of view. For example, theoretical groups 
have recently pointed out that the cooling of hot carriers 
in doped graphene is very slow when a majority of the 
carrier distribution is below the optical phonon energy. 
In this case, car rier c ooling can occur only via acoustic 
phonon emission^ 17 * 18 *. Compared to other semiconduc- 
tors, graphene stands out due to its rather large optical 



phonon energies (~ 0.196 eV and 0.162 eV). Therefore, 
the optical phonon energy bottleneck in carrier cooling 
is expected to play an important rol e in m any graphene- 
based electronic and optical devices^^. In this paper, 
we report observations of this bottleneck in the cooling 
of photoexcited carriers for the first time. 

Ultrafast terahertz spectroscopy is a useful tool to 
study the relaxation dynamics of the low energy carriers 
near the Dirac point. Previous studies of graphene using 
optical-pump terahertz-probe spectroscopy, carried out 
at room temperature by the authors^ and others^, at- 
tributed the observed relaxation transient occurring over 
a 1-10 ps time scale to carrier recombination. However, 
recent theoretical results reported by the authors and 
others show that the interband recombination and gener- 
ation mechanisms in graph ene, such as Auger scattering 
and impact ionizatio n^ 1 * 22 !, optical phonon scattering^, 
and plasmon scattering^, can have characteristic times 
much shorter than one picosecond. In particular, plas- 
mon scattering can be extremely fast, with time scales 
on the order of a few hundred femtoseconds^. In this 
paper, we present results from optical-pump terahertz- 
probe spectroscopy of photoexcited carriers in graphene 
at different temperatures. We vary the temperature in 
order to better understand the role played by several scat- 
tering processes in the observed relaxation transients. 
Our results show that the tails of the relaxation tran- 
sients, as measured by the differential probe transmis- 
sion, become remarkably slow at low temperatures. They 
extend well beyond several hundred picoseconds at tem- 
peratures below 50K. 

We also present a theoretical model to explain the mea- 
sured data. Our model includes intraband carrier scat- 
tering from optical and acoustic phonons as well as inter- 
band carrier recombination and generation from optical 
phonons, plasmons, Auger scattering, and impact ioniza- 
tion. Our model shows that the photoexcited electrons 
and holes equilibrate with each other within one picosec- 
ond due to the very fast recombination and generation 
processes. This fast equilibration causes the electron and 
hole Fermi levels to merge. The experimentally observed 
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relaxation transients beyond a few picoseconds are then 
due entirely to the cooling of the carriers. We attribute 
the very slow tails of the relaxation transients observed 
at low temperatures to the bulk of the carrier energy dis- 
tributions moving close enough to the Dirac point such 
that both intraband and interband scattering of carriers 
via optical phonon emission become inefficient. For com- 
pletely undoped samples, with symmetric electron and 
holes distributions, this occurs when the bulk of the elec- 
tron and hole distributions move below half the optical 
phonon energy. O ur res ults compliment the earlier the- 
oretical prediction d 17 ! 18 !, and our model agrees very well 
with our measurements at all temperatures. 

Graphene samples used in this work were grown epi- 
taxiall y v ia thermal decomposition of the Sl-SiC (0001) 
surface 2 ^. The samples were characterized with Raman 
and optical transmission spectroscopy to estimate the 
number of graphene layers. In experiments, the sam- 
ples were placed in a Helium cryostat. Optical pump 
excitations were performed with ~90 fs optical pulses, 
with 780 nm center wavelength, obtained from a 81 MHz 
Ti:Sapphire laser. The optical pump pulses had maxi- 
mum energies of ~11.4 nJ and were incident on the sam- 
ple from an angle with a spot size on the sample of ~1.0 
mm 2 . The maximum estimated photoexcited carrier den- 
sity was ~ 5 x 10 10 1/cm 2 . The non-equilibrium carrier 
distribution in the graphene layers was probed by moni- 
toring the differential transmission of few-cycle terahertz 
pulses, with a peak frequency of ~1 terahertz, focused to 
a 0.5 mm 2 spot size on the sample. The terahertz pulses 
were generated and detected in a terahertz time-domain 
spectrometer setup based on photoconductive antennas^ 
and had a power SNR I 10 6 . The temporal resolution of 
our pump-probe measurement was limited not by the du- 
ration of the terahertz pulses but by the duration of the 
optical pulses to ~1 ps, considering the fact that the opti- 
cal pump was incident on the sample from an angle. The 
optical pump and terahertz probe beams were chopped 
at 333 Hz and 400 Hz, respectively, and the differential 
transmission signal was detected using a lock-in amplifier 
at the sum frequency. 

The frequency dependence of the graphene conduc- 
tivity has been shown to have a Drude like behavioi!^. 
Changes in a terahertz probe pulse from the graphene 
layers are in general complex. However, in cases where 
the carrier momentum scattering time r is much shorter 
than the time scale associated with all other dynamics 
of interest, then the transmitted terahertz pulse is essen- 
tially the input terahertz pulse with an amplitude mod- 
ulation given by the instantaneous carrier distribution 19 . 
FigjT] shows the field amplitude E of a terahertz pulse 
transmitted through a 14 mono-layer (ML) graphene 
sample with and without optical pumping. Their differ- 
ence, AE, is also shown. The maximum observed phase 
shift is less than 50 fs, indicating that the electronic scat- 
tering time r is indeed very short in our samples. There- 
fore, we track only the peak amplitude of the terahertz 
probe pulse as a function of the probe delay. 




THz Sampling Delay (ps) 

FIG. 1: The electric fields of terahertz pulses transmit- 
ted through epitaxial graphene at T su b = 18K are plotted 
with and without optical pumping. Also shown is AE = 
(E pumped) — (E no pump). A small phase shift of less than 
50 fs between the two pulses indicates that the Drude mo- 
mentum scattering time, t, is small. 



The normalized differential transmission amplitude 
AT/T of the terahertz probe pulse, as a function of the 
probe delay time, is shown in Figj2]for different substrate 
temperatures for the same 14 ML graphene sample. The 
following observations can be made from this data. Im- 
mediately after photoexcitation, AT/T decreases on a 
time scale of ~1 ps, although this observation is limited 
by the resolution of our setup. After ~1 ps, AT/T recov- 
ers over time scales that strongly depend on the substrate 
temperature T su b. Decreasing the substrate temperature 
causes the recovery times to increase from tens of picosec- 
onds at room temperature to hundreds of picoseconds at 
low temperatures. At low substrate temperatures, two 
distinct time scales are observed; a fast recovery phase 
lasting to about 50 ps and a much slower phase lasting 
to hundreds of picoseconds. Measurement of the tran- 
sients over time scales longer than 300 ps were not possi- 
ble with our setup. Also, the peak magnitude \AT/T\ is 
larger at lower substrate temperatures, increasing by as 
much as an order of magnitude at T su b — 18K compared 
to T sub = 300K. 

In this section, we discuss our model, which is in good 
agreement with the measured data. We assume that im- 
mediately after photoexcitation, the photoexcited elec- 
trons and holes thermalize among themselves as well as 
with the existing carriers through carrier-carrier scatter- 
ing on a time scale of tens of femtoseconds^. After 
thermalization, the carriers have Fermi-Dirac energy dis- 
tributions with well defined Fermi- levels (Et e and Efh 
for electrons and holes respectively) and a common tem- 
perature T e h. We simulate the subsequent evolution of 
the electron, hole, and phonon distributions using rate 
equations for the electron density n e , hole density rih, 
electron/hole temperature Teh, optical phonon occupa- 
tion numbers nr and Uk at the Y and if -points of the 
Brillouin zone, and the acoustic phonon temperature T a . 
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FIG. 2: (a) Measured differential terahertz probe transmis- 
sion AT/T is plotted as a function of the probe delay for 
a ~14 layer epitaxial graphene sample at six different sub- 
strate temperatures (7^6=300, 200, 150, 100, 50, and 18K). 
The estimated photoexcited carrier density is 5 x 10 10 1/cm 2 . 
Lower substrate temperatures result in larger peak \AT/T\ 
values and slower relaxation rates, (b) The same data as in 
(a) shown on a log scale. 



The intraband and interband electron and hole scattering 
via optical phonons is described according to the mod- 
els presented by Wang et. alP^l and Rana et. alP^. The 
interband electron and hole scattering via Auger and im- 
pact ionization is described using the model of Rana et. 
al.^. The intraband electron and hole scattering via lon- 
gitudinal acoustic phonons is described according to the 
model presented by Suzuura and AndcP using a value of 
19 eV for the deformation potential. The interband elec- 
tron and hole scattering via plasmons is described using 
the model of Rana et. alP2. The anharmonic decay of 
optical phonons into acoustic phonons is described phe- 
nomenologically with the time const ant T OP t, which has 
typical values in the 0.5-2.5 ps rang d 11 | 16 | ^ J The loss 
of heat from the graphene acoustic phonons into the sub- 
strate is also described phenomenologically with the time 
constant r s „f,, which has values in the 25-200 ps ranged 

A feature of these rate equations is that the rate of 
carrier temperature change, dT e h/dt, is related to the 
rate of change of the carrier energy density, dU e /dt + 
dUh/dt, and carrier densities, dn e /dt and dn^/dt, via the 
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FIG. 3: Simulation results of the electronic temperature, T e h, 
the optical phonon temperature, Tr(« Tk), and the acoustic 
phonon temperature, T a , as a function of time after photoex- 
citation. T su b = 18K. 
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Here, C e and Ch are the electron and hole heat capacities, 
and Xe is given in terms of integrals of the Fermi-Dirac 
distribution /(.) as 



dEE 2 f(E fe ,T eh )(l- f(E fe ,T eh )) 
r °° dEEf(E fe , T eh )(l - f(E fe , T eh )) 



(2) 



The expression for Xh is obtained by substituting —Efh 
for Ef e in the above equation. Eq[l]and Eq[2]show that 
recombination (generation) due to carrier-carrier interac- 
tions (Auger scattering, impact ionization, and plasmon 
emission and absorption), in which the total energy of 
the electrons and holes does not change, always results in 
an increase (decrease) in the temperature of the carriers. 
Also note that for all the recombination and generation 
models we consider, when Efh > Ef e , the generation 
rate exceeds the re combin ation rate. When Efh < Et e , 
the opposite is tru c 1 21 * 23 * 2 -^. 

Representative results from the simulations are plotted 
in Fig(3]and Fig(4] Figj3] shows the carrier temperature 
T e h, optical phonon temperatures Tp(« Tk), and the 
acoustic phonon temperature T a plotted as a function of 
time. FigQa) shows the electron and hole Fermi lev- 
els, Ef e and Efh, and FigQb) shows the electron and 
hole densities, n e and n^. In simulations, the photoex- 
cited carrier density was assumed to be 5 x 10 10 1/cm 2 , 



18K, t = 28 fs, r, 



Opt 



1.6 ps, r sub 



100 ps, 



and the sample was assumed to be slightly n-doped with 
doping density = 1 x 10 11 1/cm 2 . Immediately after 
thermalization, the carrier temperature is large (~ 2300K 
in Figj3]) and Efh » Ef e . In the initial phase lasting 
to about 0.5-1.0 ps, the carriers lose energy via optical 
phonon emission and their distribution cools down rel- 
atively fast. In this phase, Efh exceeds Ef e , so car- 
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FIG. 4: (a) Simulation results of the Fermi levels, Ef e and 
Efh, for electrons and holes, respectively, as a function of time 
after photoexcitation. (b) Simulation results for the electron 
and hole densities, n e and nu, as a function of time, (c) 
Simulation results of the differential terahertz transmission 
\AT/T\ transient. Dashed lines indicate the approximate rate 
of relaxation at short and long time scales. 



rier generation rate exceeds carrier recombination rate 
and the electron/hole densities increase. As discussed 
earlier, carrier generation via carrier-carrier interactions 
also contributes to a decrease of the carrier temperature. 
Within a picosecond, the electron and hole populations 
equilibrate and their Fermi levels merge. Beyond ~1 ps, 
the electron and hole distributions can be approximately 
described by a common Fermi level; but strictly speak- 
ing now Ef e > Efh, so recombination exceeds generation 
and the electron/hole densities decrease with time. In the 
same time frame of 0.5-1.0 ps, optical phonon emission 
leads to an increase in the optical phonon temperatures 
until the carrier and the optical phonon temperatures 
become nearly identical. After this point, hot optical 
phonons become the main bottleneck for further carrier 
cooling and carrier cooling slows dramatically 11 , as in- 
dicated in Figj3j As the electron and hole distributions 
cool further, at some point the bulk of these distribu- 
tions moves close enough to the Dirac point such that 
both the intraband and the interband scattering of car- 
riers via optical phonon emission become inefficient in 
cooling the carriers. For the values considered in the sim- 
ulations presented in Fig|3]and Fig|4j this occurs when 
the carrier temperature falls below ^250K. Since longi- 
tudinal acoustic p honon scattering is also inefficient in 
cooling the carriers^^, the carrier cooling rate slows 
down further. This optical photon energy bottleneck is 
indicated in Fig|3j Note that the larger heat capacity of 
the acoustic phonons results in the maximum change in 
the acoustic phonon temperature being smaller than the 
maximum change in the optical phonon temperature. 

The transmission of a terahertz pulse through N 
graphene layers on a SiC substrate normalized to the 



Table. 1: Fitting Parameters 
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transmission through the SiC substrate is given by the 
expression^, 

T 1 

= (3) 

T SiC 1+ JV(7(w)V(l + «Sic) 

Here, nsic is the refractive index of the SiC substrate 
and <j(ui) is the intraband conductivity of graphene given 
bjE3 

e 2 /irh 2 r°° 

a{u) = i^—- / (f(E fe ,T eh ) + f(-E fh ,T eh ))dE 
u + i/t J q 

(4) 

FigQc) shows the calculated relative differential trans- 
mission AT/T of the terahertz probe pulse as a function 
of the probe delay. In the first ~1 ps of the simula- 
tion, the graphene conductivity increases and the tera- 
hertz transmission decreases. This increase of conductiv- 
ity is due to two factors. First, as the temperature of the 
carriers decreases, the graphene conductivity as given by 
Eq|4] increases. This results from graphene conductivity 
depending on both the total number of carriers and also 
on the carrier distribution in energy. For the same num- 
ber of carriers, the conductivity is larger if the carrier 
temperature is smaller. Second, carrier generation also 
contributes to an increase in the number of carriers and, 
therefore, an increase in the conductivity. Beyond ~ 1 
ps, the conductivity decreases and the terahertz trans- 
mission increases. This decrease in conductivity is due 
to the decrease in the carrier densities as the carriers cool 
down. The relaxation of the AT/T transient exhibits two 
distinct time scales: a first fast relaxation phase lasting 
to about 50 ps during which both intraband and inter- 
band optical phonon emission is efficient in cooling the 
carrier distributions, and the second slow phase lasting 
longer than hundreds of picoseconds during which optical 
phonon emission is inefficient in cooling the carriers. 

Figj5] shows the comparison between the theoretical 
model (solid lines) and the measurements (circles) of the 
differential terahertz transmission AT/T for three differ- 
ent substrate temperatures, T suu — 300, 150, 18K. The 
only three fitting parameters used in the simulations were 
the graphene doping density N^, the carrier momentum 
scattering time r, and the optical phonon decay time 
Topt- Table. 1 gives the fitting values used for different 
substrate temperatures. The time constant T S ubi which 
describes the loss of heat from the graphene layers into 
the substrate, was varied between 25 ps and 100 ps^l 
however, its value was found to have no significant effect 
on the simulation results for AT/T. The comparison be- 
tween the simulations and the measurements is seen to 
be very good. The model reproduces the two distinct 
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FIG. 5: Measured AT/T transients are plotted along with 
the simulation fits for different substrate temperatures. 
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FIG. 6: Measured AT/T transients for ~14 ML and ~30 
ML epitaxial graphene samples are plotted for different sub- 
strate temperatures. AT/T transients for the 14 ML and 30 
ML samples are normalized to the same peak magnitude and 
shifted along the vertical axis at each substrate temperature 
for clarity. Although the sample thicknesses differ by a factor 
of ~2, the measured relaxation time scales are nearly identi- 
cal. 

time scales observed experimentally in the transmission 
recovery transients at low substrate temperatures. The 
substrate temperature dependencies of the fitting param- 
eters r and T op t (Table. 1) are in reasonable agreement 



with the expectation d 16 * 30 ^. However, the mechanism 
responsible for the temperature dependence of the dop- 
ing density N d , although small, is unclear. 

In order to explore whether heat transfer either among 
the graphene layers or between the graphene layers and 
the substrate could be responsible for any of the features 
observed in our experiments^!, we performed measure- 
ments on epitaxially grown graphene samples with differ- 
ent numbers of graphene layers. However, we observed no 
significant changes in the observed time scales between 
the samples. Fig|6] shows the results obtained for sam- 
ples having 14 and 30 mono layers (ML) of graphene at 
different substrate temperatures. For clarity in compar- 
ison, AT/T values shown in Figj6]have been normalized 
so that the peak values for the 14 ML and 30 ML samples 
are equal. The time scales associated with the transients 
are seen to be virtually identical despite the large differ- 
ence in the number of graphene layers between the two 
samples. These observations are consistent with our the- 
oretical model where, as discussed earlier, the values of 
the time constant T su b i n the 25-100 ps range are seen to 
have minimal effect on the simulation results. 

To conclude, we have studied the relaxation dynamics 
of photoexcited carriers in epitaxial graphene at differ- 
ent substrate temperatures using optical-pump terahertz- 
probe spectroscopy. The observed differential terahertz 
transmission transients show very long relaxation times 
extending out to hundreds of picoseconds at low substrate 
temperatures. We also have presented a model that is in 
good agreement with the measurements and shows that 
carrier cooling in graphene can be very slow if the bulk 
of the carrier energy distribution falls low enough that 
both the intraband and interband optical phonon emis- 
sion processes are not possible. 
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